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Non-perturbative phenomena in four-wave mixing spectra of semiconductors are studied using 
the exact solution of a widely used phenomenological non-linear equation of motion of the exci- 
ton polarization. It is shown that Coulomb interaction, included in the nonlinearity, leads to two 
characteristic effects, which are essentially of dynamical origin, — a split of the exciton peak and 
a non-monotonous dependence of the response at the exciton frequency on the magnitude of the 
external field. Relations between the spectral features and the parameters of the system is obtained. 
It is found that the transition from perturbative to non-perturbative regimes is controlled by pa- 
rameters inversely proportional to the decay rate. It implies that the condition of low excitation 
density does not necessarily warrant applicability of the perturbational approach. 
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I. INTRODUCTION 



J2 ■ 

One of the clearest manifestations of the many-body effects in semiconductors is the the phenomenon of the non- 
linear optical response. The Pauli blocking and the Coulomb interaction between the quasiparticles lead to the 

S dependence of the polarization dynamics on its spatial distribution>ii2i&^i&£ When a semiconductor quantum well 
is excited by two successive pulses, they produce the signal in directions which are prohibited in the linear regime 
while still conserve momentum in four-wave and multi-wave mixing. The dynamical origin 8 of the formation of the 
four-wave mixing signal is naturally incorporated into the description in terms of the exciton modes characterized 
by the frequency luq and the in-plane wave vector k. From this perspective the effects of the Pauli blocking and 
. of the Coulomb interaction are clearly different. The Pauli blocking reduces locally the intensity of the field-matter 
interaction according to the magnitude of the local polarization. As a result, the excitation field with particular 
value of the wave-vector q becomes coupled to the exciton modes with different k ^ q. The Coulomb interaction 
between the excitons leads to the direct coupling between the exciton modes. In particular, when two modes with 
ki and k2 are excited, the relevant coupled modes are characterized by k^") — k.2 + nAk, where — oo < n < oo is an 
integer number and Ak = ki — k2. Although initially all energy is concentrated in the modes k^ ^ and k^ during 
the evolution the energy is redistributed among the coupled modes. In particular it leads to the formation of the 
four -wave mixing (FWM) signals, which correspond to n = 2 and n = — 1. 

It follows from this picture that the redistribution of the energy between the modes characterized by different k 
should not be the only manifestation of the mode coupling. This coupling should also lead to the modification of the 
frequencies of the exciton modes. Indeed, without the interaction and neglecting the dispersion of the exciton modes 
one has a many- fold degeneracy at the exciton resonance frequency loq. The interaction between the modes should 
lift the degeneracy and, if the excitation is sufficiently strong, make the split of the exciton peak visible in the FWM 



spectrum. 



II. THE NONLINEAR DYNAMICS OF THE EXCITON POLARIZATION 



The effect of the modification of the frequencies of the exciton modes caused by the coupling is a nonperturbative 
effect In order to study the problem we present the dynamics of the exciton polarization as the exact solution to 
the phenomenological nonlinear equation^ 

P(t; r) = -(iwo + j)P(t, r) - i(3\P(t; r)\ 2 P(t; r) 

^E{t-v)[l-\P(t;v)\yPl t ], ( " ] 

where ujq is the detuning, i.e. the difference between the frequency of the external field and the exciton frequency in the 
stationary frame, 7 is the phenomenological decay rate, E(t; r) is the envelope magnitude of the external field, P 2 at is 
the exciton saturation density and, finally, f3 = (3' — if3", with (3', (3" > 0, is a phenomenological parameter quantifying 
the interaction between the excitons. The positive real and negative imaginary parts of this parameter constitute the 
excitation induced shift (EIS) and the excitation induced decay (EID), respectively. This phenomenological equation 
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FIG. 1: The semiconductor is excited by two short pulses separated by the delay time r =t"_ — t\ 



has provided valuable interpretations for nonlinear measurements . 1 1 ' 1 2 ' 13 . The terms may be viewed as the short-time 
limit of the memory function^ which is exact to the third order in the exciting electric field. 14 The nonpertubative 
solutions have been considered by Refs. [Tolfl^ and many papers referenced in Ref. [13 

We consider the excitation of the semiconductor by two short pulses (see Fig. Q]) acting on the system during the 
intervals t'_ < t < t' + and t"_ < t < t'{ , respectively. Between the pulses, t' + < t < t"_ and after the second pulse 
t > i" the dynamics of the exciton polarization is free and satisfies Eq. fl} with E — 0. We study the free dynamics 
polarization Po(t) before calculating the polarization P\^{i) of the driven dynamics. Then we are going to solve for 
the polarization in the following order: P\(t) between t'_ < t < t' + , Pio(i) between t' + < t < t'!_, ^(i) between 
t'!_ <t < t'l_, and at the end the final solution P2o(t) for t > t'{_. 

The free dynamics [E(t) = 0] can be solved exactly by noting that Eq. (JTJ) in this regime is reduced to the simple 
form 

P (t; r) = - [^o + 7 + *n(t; r)] P Q {t; r), (2) 
where Q(t; r) = 0\P o (t; r)| 2 , and by observing that Eq. ([T]) yields an equation for the magnitude of the polarization 

||P (i;r)| 2 = -2[ 7 + /?"|Po(i;r)| 2 ]|P (t;r)| 2 (3) 

with the solution 

|P (t;r)| = |P (0;r)| e -^(*;|P (0,r)|), (4) 

where 

-l 

(5) 

The free propagator may be expressed in terms of the amplitude modulation A(t; r) and the phase modulation tp(t; r), 
P o (t;r)/P o (0;r) =: 9 (t; |P o (0,r)|) = A(t; |P (0, r)\)e^ ot-jt+i^r) ^ (6) 

where 

^;r) = |^lnA(t;|Po(0,r)|). (7) 

The amplitude term shows the nonpertubative effect of EID. It is interesting to note that EID does not lead to a 
mere modification of the decay rate, 7. Instead the amplitude modulation A(t; r) decreases to a fraction of the initial 
amplitude at twice the linear rate, i.e., 27. The phase modulation shows an oscillation dependent on the nonlinear 
quality factor /?'//?". This is related to the Goldstone mode in the excitons studied in Ref. YLn 



A 2 (i;|P (0,r)|) = 



i+ ^|Po(Q;r)| 2 2 

7 
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FIG. 2: (a) Four-wave mixing spectrum corresponding to 2ki — k2 as a function of the normalized pulse area. The parameters 
of the system are chosen to be 13'/"/ = 5 ■ 10 2 , /3'/2/3" = 5. The horizontal and the vertical lines show the positions of the 
sections of the surface P^(ui; e) presented in (b) and (c). (b) The FWM spectrum at the fixed value e/Pg at = 0.19. (c) The 
dependence of the FWM response at the exciton frequency P^(u>o; e) as a function of the pulse area. 



The initial conditions for Eq. © are determined by the polarization distribution right after the external field is 
switched off. We find the immediate response of the system assuming that the duration of the excitation pulses is 
much shorter than the typical dynamical time scales determined by detuning and the decay rate. 

We turn now to the driven (E(t) ^ 0) time evolution of the polarization, i.e. we consider the time interval 
t_ < t < t + where the particular excitation pulse does not vanish and factor out the term exp(ik • r) so that 
we can consider the excitation pulse to be spatially homogeneous. Neglecting the contribution to the phase ~ 
f* + [iojQ + 7 + iQ(t; r)] Pi^(t] r)dt we solve the dynamical equation and find the relation between the polarization at 
the instances t_ and t+ 



Pi, 2 {t+) = P' h2 (t- 



p2 

sat 



P'i j2 (i_)tanh 



p2 

sat 



p'Ut-) 



p2 

sat 



atanh 



P"l,2(t-) 



p2 

sat 



P'Ut-) 



(8) 



where P'l^t-) = Re[Pi.2(£-)] and P"i.2(^-) = Ini[Pi. 2 (£-)], and we have introduced the area of the exciting pulse 
e = J dtE(t). 

Using the solutions for the free and driven polarization in Eqs. (|2|) and ([8]), we solve the time evolution of the 
polarization in all the four regions t'_ < t < t' + , t' + < t < t'!_, t"_ < t < t'L, and t > t'{_. Assuming that the system 
initially is in the ground state P(t'_) = we find the exciton polarization created by the first pulse 



where 



P 1 (t' + -v) = -ie^ r P 1 (t' + ) 



Pi(f + ) = Psat tanh 



(9) 



(10) 



It follows from this equation that the saturation effect for the first pulse reduces to a simple (although nonlinear) 
renormalization of the signal area. Using Eq. © in Eq. (01 we can see that the effective frequency f2(f; r) determining 
the dynamics of the polarization is constant across the sample. As a result there is no coupling between the exciton 
modes characterized by different k. 

Denoting the delay time, the time separation between the pulses, by r = t"_ — t' + we obtain the polarization right 
before the arrival of the second pulse 



Pio(i";r) =Pi(*' + ;r)0 lo (T), 



(11) 
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where #io(t) = 6(t; \Pi(t',)\). The spatial distribution of the polarization Pio(t";r) plays the role of the initial 
condition for the immediate response with respect to the second pulse according to Eq. ([5]). For the analysis of 
the time evolution of the polarization during the second pulse it is convenient to factor exp(ik 2 • r) out of P2(t;r) 
introducing 

P 2 (t; Ak-r) = e-* k2 ' r P 2 (i;r) (12) 

with Ak = ki — k 2 . The reduced distribution P 2 (t; Ak • r) satisfies Eq. (fTJ) with the modified spatial profile of the 
external field E(t) = e~ lk2 ' r E(t; r). The form of the initial conditions in Eq. ([8]) changes according to 

P 10 (^;Ak-r) = -i0 10 (r)P L (t' + )e iA1 '- r . (13) 

It follows from Eqs. (JSj> and (|13|) that for the second pulse the role of the saturation effect is two-fold. It modifies 
the pulse area and excites all modes k^ rather than just a single mode as we had for the first pulse. We present the 
polarization as a superposition of the multi-wave mixing modes 

P 2 (t; K ) = Y / Pt ) (ty nK , (14) 

n 

where k = Ak • r and P 2 n \t) are the amplitudes of the multi-wave mixing polarizations. Substituting this represen- 
tation into Eq. we can see that in terms of the multi-wave mixing modes the free dynamics of the polarization can 
be presented as the dynamics of a system with an infinite number of degrees of freedom coupled to each other due to 
the nonlinearity 

^P 2 W(t) = -N + 7)^(1) - i^O re _ ro (t)P 2 (m) (t), (15) 

m 

where f2 n (f) = (3(2tt)~ 1 J_ \P 2 {t; K)\ 2 e~ lK,n dn. The initial conditions for Eqs. (fT5|) are constituted by the immediate 
response to the second pulse thanks to the saturation effect. This picture clearly illustrates the difference between 
the effect of the Pauli blocking and the Coulomb interaction on formation of the multi-wave mixing response. 
Using the solution of the equation of motion for free polarization dynamics we find 



p(«) 



(t) = ^f Mt'U k)0 (t; |iMt£;/s)|) e- iKn dK. (16) 



This equation and Eqs. (|5J| and (| 1 3[) give the exact evolution of the exciton polarization in the limit of short excitation 
pulses in the two-pulses scheme. Formally one can obtain from Eq. (|16p the spectrum P^ n '(u>) using the Fourier 
transform of 8(t; |P 2 (i" ; k)\) with respect to time 

9(u; |P 2 (t + ; k)\) = i{w + i/2) 2*1 \l; 2 +^;-- iw; j , (17) 

where w = (lo — cjq)/27, X — (3'/2/3" and 2 Fi is the hypergeometric function. Technically, however, because the 
second argument of this function is a complex number it may be more efficient to calculate the spectrum using the 
time series P 2 (t). 

The spectrum corresponding to the four-wave mixing directions 2k 2 — ki and 2ki — k 2 is obtained from Eq. (|16p 
taking n = — 1 and n — 2, respectively, and for n — 2 is shown in Fig. [5] as a function of the pulse area. In order 
to estimate the relation between the decay rate and the non-linear parameter we have used the expressions following 
from the microscopic consideration 7 (3' w 1.52a^£';, and P^ at = 7/47ra^, where is the exciton Rydberg and is 
the exciton Bohr radius. 

The spectrum has two typical features. The first one consists of branches detached from the exciton frequency with 
increasing pulse area. These branches for sufficiently high amplitudes of the excitation field may manifest themselves 
on the spectrum in the form of multiple resonances. The second interesting feature is the oscillatory character of the 
field dependence of \P 2 2 ^ (loq; e)| 2 , that is the FWM response at the exciton frequency. 
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FIG. 3: The form of the four-wave mixing spectrum in low saturation regime in the case of negligible EID (J3" = 0) and 
/3'/t = 2.5 ■ 10 2 . The solid line shows the position of the resonant frequency determined by uj — ljq + (31 + r/j. 



A. The FWM spectrum in the limit of negligible EID 

We start the discussion of these features from the simplest case (3" — (vanishing EID) and high P sat (low saturation 
regime). The first assumption simplifies the effect of the initial conditions on the polarization dynamics while the 
second simplifies the relation between the excitation field and the polarization of the immediate response, so that 

P 2 (t'U k) = -i (|eiMT)|e- lK+ * + |e 2 |e**») , (18) 

where 4>\ = arg(ei^i (r)) and 4*2 = arg(e 2 ). 

It is interesting to note that according to Eqs. (|15p and Q in this approximation since only Q±i(t) differ from 
zero the dynamics of the polarization is described by the equation of motion for a ID tight-binding model with the 
time-dependent coupling between neighboring sites oc /?|e2ei#io( T )|e~ 27t . The relation between the amplitudes of 
the excitation at different sites gives the relation between the amplitudes of the signals corresponding to multi-wave 
mixing. Initially the excitation is localized on sites n = 0, 1 and with time it propagates along the chain giving 

P a W (t) = -e^-^-^+^expj-iwoi-Ti-^ (l -e~ 27t ) j x ^ 

{|ei0io(r)| J n -i [t?(1 - e~ 2 ^)] + i\e 2 \J n - e^ 4 )] } , 

where J n are the Bessel functions of the first kind, I = |ei0i O (r)| 2 + |e2| 2 , and rj = P\e2eiOio(r)\/j. Deriving Eq. (fl9|) 
we have used the Jacobi-Anger expansion^ exp(zzcosK) = X)^L-oo i n Jn(z)e mK . 

The spectrum of the four-wave mixing signal corresponding to n = 2 is shown in Fig.[3J Shortly after the excitation, 
for t <C I/27, the amplitude of the multi-wave mixing signal drops exponentially with the order of mixing, ~ (j]jt/2) n . 
The exponential drop holds asymptotically in time if 77 < 1. This result agrees with the perturbational approach. The 
situation, however, drastically changes if the parameters of the system are such that rj 3> 1. In this case, the intensity of 
the multi-wave mixing signal becomes independent of its order starting time t > — ln(l — i]~ 1 )/2-f w l/2/3| 6261^10 (t)|. 
This consideration suggests naturally to identify 77 = 1 as a critical value that separates perturbative and non- 
perturbative regimes. 

The detailed form of the spectrum is essentially determined by the fact that the effective coupling between the 
modes vanishes with time. Right after arrival of the second pulse the polarization oscillates with the frequency 
detuned from the exciton frequency by the value ~ (31 + 777, where the second term estimates the contribution of the 
Bessel functions. This detuned frequency qualitatively describes the dependence of the frequencies of the detached 
resonances on the signal area (see the bold solid line in Fig. [3]) . In particular it shows that the resonances in the 
spectrum widen with increasing nonlinear parameter (3. In the opposite limit, t ^> 1/7, the coupling between the 
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modes vanishes and, as a result, the exciton polarization oscillates with the non-modified exciton frequency uiq. These 
oscillations give rise to the resonant behavior at the exciton frequency, for which dependence on the pulse area is 
determined by the asymptotic value of the Bessel functions J{rj). Therefore, the response at the exciton frequency 
essentially depends on whether the system is in perturbative or non-perturbative regime. We illustrate the difference 
between these regimes considering the effect of the positive and negative delay time. The expression for the negative 
delay time can be obtained from Eq. (|19|) by considering the signal in the "conjugate" direction n = 1 — n and 
exchanging e% <-> 62- 

In the perturbative regime, ij C 1, we obtain 



P 2 (n) (o;o;r < 0) ~ Arf^ 2 Wnfffc 1 {-t)\ 2 , 



(20) 



where rjo — /5|eie2 I n this regime l^ioMI ~ ex P( — 7 r )j thus, the positive delay signal decays with time constant 
2wy while for the negative delay it decays more slowly with the constant 2(n — 1)7. For the case of FWM signal this 
result agrees with the perturbational calculations^ 

In the non-perturbative regime the positive and negative delay signals are determined by the oscillating asymptotics 
of the Bessel functions^ 



1 

nr] 



(-1)" (|eiMr)| 



C2| 
2 



£2 



sin(2?7) 



(21) 



where "— " ("+") sign corresponds to the positive (negative) delay. Writing Eq. (f2Tj) we have neglected the oscillating 
term oc cos(2?7) vanishing as (|ei0io(r)| + [62] )/rj 2 with increasing signal area. Thus, for both positive and negative 
delays the response at the exciton frequency saturates at the oscillations with the period r\T — 27r. The strong 
asymmetry between these cases specific for the perturbative regime does not hold any longer and the only difference 
is the phase of the oscillations. It should be noted that in Ref. [2(J the saturation of the FWM response was attributed 
to the renormalization of the pulse area by the EID and EIS. The present consideration, however, suggests that the 
origin of the non-trivial dependence of the response on the pulse area is the redistribution of the excitation over 
multi-wave mixing modes. Considering the identity J2 n J 2 {z) = 1 for the limits z«l and z 3> 1 it can be seen that 
such redistribution is especially effective in the nonperturbative regime resulting in essential suppression of the FWM 
response. 

It should be emphasized that the mechanism of the oscillatory dependence of the response at the exciton frequency is 
different from the Rabi oscillationS ) 21 ' 22 which would correspond to the non-monotonous dependence of the immediate 
response on the excitation field. In the case under consideration the oscillations are the result of the free dynamics 
of the exciton polarization when the external field is turned off. The physics of the Rabi oscillations and of the 
oscillations of |P 2 (u; rj)\ 2 , of course, are essentially the same. As noted above the dynamics of the polarization in 
the case under consideration appears analogous to a ID tight-binding model with vanishing coupling between the 
neighboring sites. In quantum mechanical terms it can be described as a multiple level system, where the levels 
correspond to the exciton modes, with time dependent field Vy(i), which couples different levels. Depending on the 
"area" of the off-diagonal elements, J dtVij(t), one has the oscillations of the final populations of the different levels. 



Translated to the language of the multi-mixing signals P 2 this result implies the oscillations of |P 2 ™ ; (cJo)| 2 since 
asymptotically, as has been noted, one has the dynamics determined by the non-perturbed exciton frequency. 



)(»)/ 



B. The effect of the excitation induced decay 



In order to study the effect of the EID on the spectrum (compare Figs. and [3]) we use Eq. (fT7|) assuming that 
|ei^io ("7") I = = e - Considering the asymptotic values of the hypergeometric function 2 ^, in the limit 77" = 2/3"e 2 /j 
1 we can approximately present the spectrum in the form 



p 2 <*>M 



C{u-X)A 2 e~ lXl ^" + — 



r)" iX -1/2 



(22) 



where 



(23) 
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and A\ = A n (l), A2 = A n (l/2) with j4 n (p) depending on X only 



A»(p) 



1 

2^ 



-ZK \ n -JTiK 



d K (l + e- lK )e 
x exp I — (iX + p) In 



2 cos' 



(!)]}■ 



(24) 



Similarly to the case of negligible EID, the response oscillates and reaches the saturation in the high excitation 
limit 77" 3> 1. At the exciton frequency the magnitude of the signal is 



^|A 2 | 2 tanln%Y). 



(25) 



This saturation also is of the dynamical origin since we work in the regime of weak saturation due to the Pauli 
blocking (e/-P 2 at < 1). The significant difference with the previously considered situation is that now the minima of 
the spectrum as a function of the pulse area are not equally spaced as the oscillating part has the form oc cos(Y In 77"). 
Thus, with the decreasing ratio between the real and imaginary parts of the non-linear parameter the crossover from 
linear to logarithmic scale occurs. 

Finally, we discuss the effect of the saturation parameter. Qualitatively, this effect can be understood as follows. The 
non-monotonous behavior of the semiconductor response on the pulse area studied above is supported by unrestricted 
increase of the polarization of the immediate response [see Eq. (|18|) ]. However, the saturation effect renormalizes the 
pulse area so that the magnitude of the polarization can not exceed P sa t- For example, from the perspective of the 
discussion of the effect of EID this means that Eq. (j2"2"|) remains valid only if two restrictions are met 2/3"e 2 /7 ^> 1 



and e < P sa t- This imposes the restriction of the decay rate to be sufficiently small 7 <C (3"P 2 at . 



III. CONCLUSION 



We have studied non-perturbative effects in four-wave mixing spectra of semiconductors. These effects are analyzed 
using the exact solution of the non-linear equation of motion of the exciton polarization taking into account excitation 
induced shift (EIS), excitation induced decay (EID) and the saturation effect phenomenologically. We found that the 
interaction between the excitons accounted by EIS leads to two specific spectral features — a split of the exciton peak 
and a non-monotonous dependence of the response at the exciton frequency ujq on the magnitude of the external field. 
The important characteristic of the splitting is that new spectral features should appear at frequencies higher than 
luq. This allows one to make a distinction between the effect of interaction of the exciton modes and the manifestation 
of bound biexciton states, which should modify the spectrum at frequencies lower than luq. 

We would like to emphasize that these effects do not appear in any order of the perturbational (x ) approach. It 
can be shown that the appearance of additional spectral features can be traced as a divergence of the perturbational 
series. It should be stressed out that the crossover from the perturbative to non-perturbative regimes is governed 
by parameters that essentially depend on the decay rate (e.g. /3|ei^io(r)e2|/7 ~ 1 in the case (5" — 0). This 
means that the condition of low excitation itself does not necessarily warrant the validity of the perturbation theory. 
As an ultimate example one can consider the model with 7 = when the spectrum (for r = 0) has the form 
P(w) oc \JyJ (uj — ujq — f3I) 2 — (3 2 I 2 with the exciton peak being splitted for arbitrary low excitations. 
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